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Effective boundary conditions for dense granular flows 
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We derive an effective boundary condition for granular flow taking into account the effect of the 
heterogeneity of the force network on sliding friction dynamics. This yields an intermediate boundary 
condition which lies in the limit between no-slip and Coulomb friction; two simple functions relating 
wall stress, velocity, and velocity variance are found from numerical simulations. Moreover, we show 
that this effective boundary condition corresponds to Navier slip condition when GDR MiDi's model 
is assumed to be valid, and that the slip length depends on the length scale that characterises the 
system, viz the particle diameter. 
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Granular media exhibit a wide range of flow regimes 
[l[ , as well as a plethora of dynamical instabilities [If. Fo- 
cusing on gravity (or shear) driven flows, three regimes 
have been pointed out: (1) the collisional {gas-like) 
regime, where energy is dissipated by the inelasticity of 
the collisions, (2) the dense flowing {liquid-like)vegime, in 
which particles undergo long lasting contacts, and dissi- 
pation occurs through dynamic friction, and (3) the static 
{solid-like) regime, which is capable to mantain struc- 
tures due to the threshold, non-linear nature of static 
friction. These regimes were studied with both experi- 
ments and discrete models, the latter having experienced 
a great advance in the last years, starting from the work 
of Cundall and StrackQ. 

Reliable continuum models would be of great advantage 
in simulating granular media, particularly when dealing 
with complex geometries or flows; in fact a unifying the- 
ory is still lacking. In this perspective, regimes (1) and 
(3) have been worked out with some success in a vari- 
ety of theoretical studies, respectively with the kinetic 
theory of granular gases [J] and with continuum critical 
state soil mechanics [![. For the dense regime, various 
theoretical approaches have been developed (and exten- 
sively reviewed in Q); the last, more attractive one, is 
that proposed by the GDR MiDi based on the inertial 
number I (see Q Q [I] [3 ) , the importance of which was 
already stated by GoddardjllJ. However, despite the 
great attention towards continuum models and Theolo- 
gies, little effort has been devoted to develop realistic 
boundary conditions for the velocity field at smooth or 
rough walls, even if the crucial role of side walls was rec- 
ognized, for example, for inclined chute flows [p|. A 
common experimental approach developed to overcome 
this issue is the practice of gluing particles to the walls, 
in order to assume a no-slip boundary condition in the 
interpretation of the results. This intelligent choice is of 
fundamental importance but has, in our opinion, two ma- 
jor drawbacks: at first, it is knownQ that for high shear 
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rates particles undergo strong slip at the glued particles 
- bulk particles interface, a slip that adds some difficulty 
in holding the continuum hypothesis; thus it is not clear 
whether the glued particles are part of the bulk or of 
a bumpy wall, so that boundary conditions must be ex- 
pressed on the first moving layer in contact with the glued 
one. The second drawback of this experimental practice 
is the partial applicability to real situations: the flow 
on smooth surfaces such as in hopper discharge usually 
shows particles slipping at the solid interface. Slip can 
be promoted or can be an undesired phenomenon, often 
we are concerned with stick-slip phenomena fl3j ■ which 
are common in dry-friction dynamics[14] ; in all of these 
cases, a deeper understanding of the behaviour of gran- 
ular materials flowing near a boundary is needed, and 
the no-slip boundary condition is not the most valid ap- 
proach. 

In a recent work[l5| we used the mixing-length model 
proposed by the GDR MiDi and showed that using a 
slip boundary condition instead of a no-slip one consid- 
erably improved the predictions of the model in the verti- 
cal chute configuration; there we used a Coulomb friction 
condition, which could be a valid alternative to the no- 
slip condition. In this work we go even further, showing 
for a simple case that taking into account the effect of 
the heterogeneity of the force network yields an interme- 
diate boundary condition which lies in the limit between 
no-slip and Coulomb friction; moreover, we show that 
this effective boundary condition corresponds to Navier 
slip-length condition if GDR MiDi's model is assumed, 
and that the slip length depends on the length scale that 
characterises the system, viz the particle diameter. We 
consider a single particle of mass m and diameter d ly- 
ing on a plane, moving with instantaneous velocity V; 
the particle is subjected to a normal force P and to a 
tangential force T. We will neglect, for simplicity, the 
effect of couples acting on the particle considering only 
traslational, sliding movements. We will assume that due 
to the heterogeneous nature of the medium the normal 
force P is a random function of time with a given distri- 
bution function. Even T could fluctuate, but we assume 
for the sake of simplicity that only the normal force does; 
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qualitative results are not affected by the choice of the 
fluctuating force. Let F be the friction force; we consider 
the simplest model of solid friction, e.g Amontons law, 
with only one friction coefficient /i: 



F = 



T if V = 0andT < [iP 
jiP else 



(1) 



The motion of the particle is calculated from Newton's 




FIG. 1: Schematism of the variables considered in this work. 



law: 



dV T „ 

m—— = 1 - b 
dt 



(2) 



If the normal force was constant, only two situations 
would be possible, corresponding respectively to no-slip 
and Coulomb conditions. But if the force fluctuated, 
the particle would undergo slip and no-slip events, which 
globally represent a non-Coulomb slip phenomenon; our 
aim is to derive an average expression for the slip veloc- 
ity as a function of the forcings. Let's consider a typical 
distribution of normal forces of the form: 



p(f) = a(l-be-' )e 



(3) 



as suggested in (l6j. This distribution of forces holds 
for normal forces in uniaxial compression, in a spatial 
sense; we will make the hypothesis that this distribution 
acts also between successive rearrangements of tangen- 
tial forces in time. Our choice is supported by the fact 
that results do not depend on the particular choice of 
the distribution, apart from one point (the existence of 
a cutoff value in the force) which will be discussed later, 
and whose influence is limited. We suppose further that 
the force is a piecewise linear function whose nodes are 
extracted from this distribution. Let P aV e be the av- 
erage value of the normal force. We choose the time 
step between successive force rearrangements to be equal 

to the relaxation time r = y' pr-^-Q; it follows directly 

that rescaling t by t the time step over which the force 
rearranges is 1. Further rescaling leads to the dimen- 

sionless variables: V = V 'J ' p m d and T' = T//j,P ave , 
P' = P I Pave- If we define a(t) as: 



Oif V 
1 else 



OandT < fiP 



the equation of motion becomes 
dV 

IF 



a t i (T' - P') 



(4) 



(5) 



from which we can compute the average rescaled slip ve- 
locity defined as: 



Ce = M lim - / dt' [ a (T' - P') dt 
^+°° T Jo Jo 



(6) 



We solve numerically the equation of motion; an example 
of the stick-slip behavior of the system is given in figure 
[2] An initially motionless particle can start to move 
only if the instantaneous normal force is below the yield 
threshold. A moving particle can decelerate only if the 
normal force is higher than the threshold. Moreover, 
it is clear from Fig. [2] that the area in which normal 
forces oppose motion is larger than the area in which 
they promote motion; it is the dynamical nature of the 
system that causes the body to have a non null average 
velocity. It would be desirable to find a relationship 
between the average slip velocity computed by means of 
Eq. [6] and the rescaled average tangential force (which 
corresponds to a rescaled effective friction coefficient, 
being fi e ff/fi — T/[iP ave ). After solving Eq. [3] it is 
possible to look at the dependence of the statistics of 
the particle motion on the average value of the force 
in Fig. [3] The curves evidence a no-slip limit at low 
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FIG. 2: Example of the local dynamics of the system. 

values of the rescaled force T', and a Coulomb limit 
for T' — > 1. The way V approaches depends on the 
nature of the distribution p(T): if the distribution had 
an upper cut-off value, it would be easy to conclude that 
the system had a sort of yield-stress behavior at the wall, 
with a finite range of T" giving V' = 0; in the other case, 
without cut-off, the average velocity would be only for 
T" = 0. This is the only point in which the choice of the 
distribution function qualitatively changes something in 
the results; however, the fast decrease of the tail in the 
distribution, if not giving a "plastic" behavior, would 
give some sort pseudo-plastic behavior, because of the 
need to impose a certain stress to obtain an appreciable 
slip. So, with a certain loss of exactness, it is possible to 
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FIG. 3: Dependence of statistics of particle velocity on statis- 
tics of force, (top) Rescaled average slip velocity vs average 
pulling force, (bottom) Rescaled velocity variance vs average 
slip velocity. Best fits from equations 1718 1 are also included. 



assume also a yield-stress formulation for the BC. 
It is interesting to note that also the variance of the 
distribution of the instantaneous particle velocities, that 
corresponds to the concept of granular temperature, 
which we express as 9 =< (V(t) — V ave ) 2 >, grows 
when T' increases; so granular temperature at the wall 
has a behaviour similar to that of the slip velocity. 
Due to its definition, 9 is made dimensionless with the 
position 6' = Ojrj- ^ n fig ure IS! correlation between 
granular temperature and velocity is shown to follow 
a power-law behavior. From a general standpoint, the 
boundary conditions can be expressed with the help of 
the following fitting functions: 



T' = 



V 



V 



Ci 



9' = c 2 V 



'0 



(7) 



(8) 



where < a, f3 < 0.5 and Ci,C2 > are fitting parame- 
ters. A very good fit is obtained for a « 0.28, (3 « 1.5, 



Ci « 0.51, C2 =~ 1.8. The fit is the solid line in the 
figures. Eqq. [7]and[5]are the simplest expressions for the 
effective boundary conditions that can be applied at the 
wall characterised by a particle-wall friction coefficent \i. 
Navier boundary condition, relating the slip velocity and 
the gradient of the velocity normal to the boundary via a 
slip length A, is a common way to characterise slip in fluid 
flows in micro and nanochannels; however, there is not 
a single plot of this condition in the V vs T' diagram, 
because such a plot needs information on the relation- 
ship between stresses and deformation rates in form of 
constitutive relations. For a newtonian fluid, 



V' = X 



mP 



r 



(9) 



which is linear and parametric in ^ w . A Bingham 
yield-stress fluid will have an explicit relation of the form: 



V' = X- 



mP 



(T' 



{for 



T' > Ty) 



(10) 



where Ty is the rescaled yield stress and rf is the viscosity 
coefficient in Bingham's model. So a Navier condition for 
Bingham's model in the V vs T" plot is a line shifted by 

Ty and again parametric in ^ \J~^f~- 
Both of these relationship obviously do not conform to 
the behavior obtained from the model developed in this 
work; assuming a mixing length model as GDR MiDi's, 



where -p — /i(/), with I = 



(the difference in the 



y/P/rnd 

expression of I of the previous literature is due to the 
fact that here P is a normal force, not a pressure), the 
assumption of fj,(I) = fi s + (taken from Pouliquen 

and coworkers [12, @|) yields for a Navier BC: 



V 



X T' - fj, s //j, 



{for r>ii./n) (11) 



which reaches an asymptote for T' — > ^.2/ ^, and is 
in the range — A* s /m- Thus, to unify the curves and 
represent the results obtained from the simple model of 
wall friction presented in this Letter, A must be a function 
of the form: 



A = k d C(T') 



(12) 



where C(T") accounts for the change in the position of 
the asymptote and can be expressed simply as: 



c 



- T' 
1 -T> 



(13) 



An important result is given in Eq. 1121 to unify the 
curves as obtained in the "experiments" , A must be a 
multiple of d: this is actually an important result, being 
d the only internal length scale of the system, and so the 
best choice as a basis for estimating the slip length. The 
typical form of V' vs T' curves for the various models is 
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FIG. 4: Rescaled average velocity vs average pulling force for 
different BCs/constitutive laws. The slope of Newton and 

Bingham lines is A-^-i/ 

given in figure [4j 

To resume, the intermediate, efficient boundary condi- 
tion we are looking for can be qualitatively expressed as 
Navier slip condition in a mixing length framework, the 
slip length corresponding to a multiple of the particle di- 
ameter. A step further can be made in the direction of 
determining a value for A. Let's admit the yield-stress 
behavior of Pouliquen's form for and suppose that 

fi2 ~ A* (remember that /1 is the particle-wall friction 
coefficient). In this perspective the slip length is sim- 
ply proportional to the particle diameter and the best fit 
gives X/d « 0.2. This value gives a sort of minimum slip 
length; in the case with [12 > fi the slip length diverges 
for T' — > 1. 

It is important to note that V' is an analog of the iner- 
tial number / for near wall flows, and T" is an effective 
wall friction coefficient as is for the bulk; thus it is 



interesting to note that the shape of the curve T'(V') is 
very close to that of /i(i); this can lead to some ideas on 
the origin of the effective friction coefficient in the bulk 
remembering that the effective wall friction coefficient 
derives from the assumption of heterogeneous forces. 
In this work we do not aim to define the correct functional 
form for these BCs (even if a very good fit was obtained 
for this simple case), but we want to underline that real 
boundary conditions (even in simplified setups) are not 
no-slip or Coulomb-like, and assuming one of these lim- 
iting BCs can introduce errors in the physical validity of 
granular flow models; this slip behavior can be captured 
by a modified Navier condition, where the slip length is 
proportional to the particle diameter. 
To resume, a simple model of a particle sliding with the 
simplest frictional law on a plane has been developed in 
this Letter to determine effective boundary conditions 
for granular flows. To account for the heterogeneity of 
the medium, the particle is subjected to a random nor- 
mal force, while a constant tangential force is assumed 
for simplicity. The dynamics consists of stick-slip events, 
which are related to the heterogeneity in the stress field; 
we reported the resulting dependence of the average tan- 
gential force on the average slip velocity and on the vari- 
ance of the velocity of the particle (i.e. granular tem- 
perature), thus providing two possible effective bound- 
ary conditions for the velocity and granular temperature 
fields. The results are well fitted by simple laws and rep- 
resent for the velocity field an intermediate behavior be- 
tween Coulomb's law (at high velocities) and the no-slip 
boundary condition. Granular temperature is related to 
the velocity by a simple power law behavior. In addition, 
we demonstrated that the curve obtained by numerical 
simulation satisfies a modified slip-length Navier bound- 
ary condition within a mixing length model of granular 
flow, with the slip length being proportional to the char- 
acteristic length of the system, the particle diameter. 
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